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Abstract 

Micro-Pattern Gaseous Detectors (MPGDs) such as Micromegas or GEM 
are used in particle physics experiments for their capabilities in particle 
tracking at high rates. Their excellent position resolutions are well known 
but their energy characteristics have been less studied. The energy resolution 
is mainly affected by the ionisation processes and detector gain fluctuations. 
This paper presents a method to separetely measure those two contributions 
to the energy resolution of a Micromegas detector. The method relies on 
the injection of a controlled number of electrons. The Micromegas has a 
1.6-mm drift zone and a 160-//m amplification gap. It is operated in Ne 
95 %-iC4Hio 5 % at atmospheric pressure. The electrons are generated by 
non-linear photoelectric emission issued from the photons of a pulsed 337-nm 
wavelength laser coupled to a focusing system. The single electron response 
has been measured at different gains (3.7 10 4 , 5.0 10 4 and 7.0 10 4 ) and is 
fitted with a good agreement by a Polya distribution. From those fits, a 
relative gain variance of 0.31±0.02 is deduced. The setup has also been 
characterised at several voltages by fitting the energy resolution measured 
as a function of the number of primary electrons, ranging from 5 up to 210. 
A maximum value of the Fano factor (0.37) has been estimated for a 5.9 keV 
X-rays interacting in the Ne 95 %-iC4Hio 5 % gas mixture. 
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1. Introduction 



Thanks to their very good resolution in position and capacity to work at 
high counting rates, Micro Pattern Gaseous Detectors (MPGD), like GEM 
[it] and Micromegas 0], have been studied to replace wire chambers for 
tracking in high energy physics experiments 0, 0> H] • More recently, those 
detectors have been used in nuclear physics experiments with heavy ions at 
energies ranging from several hundred keV up to a few MeV per nucleon. For 
example, a TPC (Time Projection Chamber) equipped with a four-GEMs 
structure was constructed for the study of the two-proton radioactivity [g] . 
For future experiments with nuclear radioactive beams, R&D activities have 
started on a TPC where the gas is both the target and the detection medium. 
This kind of detector must fulfil some specific conditions like the possibility 
of varying the pressure from a few tens of millibars to a few bars. An 
important feature is also the energy resolution that could be obtained on a 
large range of primary electrons depending on the ionisation power of the 
nuclei involved in the reactions. A small prototype with an active area of 
16 cm 2 has been constructed for dedicated studies on position and energy 
resolutions in the specific field of detection in nuclear physics. It is equipped 
with a Micromegas of a fairly large amplification gap (160 //m) in order to 
reach a sufficient gain at low pressures [tJ. Before starting the studies at 
different pressures with different gas mixtures, we decided to investigate the 
energy resolution in more standard conditions, i.e. at atmospheric pressure 
and with a gas mixture that could allow a high dynamic range of gains 
without sparking. Different tests with Argon based gas mixtures showed that 
it is not possible to reach a sufficiently high gain with our large amplification 
gap and we chose Neon with a low amount of quencher: Ne 95 % iC4Hio5 %. 
The energy resolution depends on the fluctuations of the number of primary 
electron-ion pairs and the avalanche multiplication in the gas. The ionisation 
fluctuations can be quantified by the Fano factor which depends on the 
gas and the energy of the primary particle ([l3]-[la]). The gain fluctuations 
depend on the gas, the detector geometry and the applied electrical field [171 . 



la]. A powerful method to measure the gain fluctuations relies on knowledge 



of the detector Single Electron Response (SER). Then, by measuring the 
charge resolution for different numbers of primary electrons, the fluctuations 
of the primary charge can be estimated and a Fano factor deduced. The 
experimental set up is described in section [2j The method for the gain 
calibration independent of the SER measurements is presented in section [3l 
It relies on the laser and the X-rays from a 55 Fe source. In section dj the 
SER at three different gains are shown together with their associated fits by 
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Polya distributions. The factors affecting the energy resolution measured 
with the laser and the 55 Fe source are discussed in section [5) 

2. Experimental setup 

The experimental setup is similar to the one described in ref. and 
is shown in Fig. (TJ The test-bench is composed of a pulsed laser, an op- 
tical line, and a 3D positioning system. The laser is a 337-nm wavelength 



pulsed Spectra-Physics VSL337 |20j] and delivers pulses with a 4-ns FWHM 
duration, an energy of 120 /iJ, a pulse to pulse energy stability better than 
4 % in standard deviation, and a repetitive rate of 10 Hz. The laser beam 
is collimated by a diaphragm and a telescope, made of two spherical lenses 
and focused by a 60-mm focal length triplet. The spot size is smaller than 
100 p. In addition, a periscope is placed between the telescope and the 
triplet for mechanical reasons to adjust the beam height. The light intensity 
is tuned by means of commercial calibrated neutral density filters. The po- 
sition and the focusing of the beam can be changed by moving the detector 
in the three directions of space with 1-fim accuracy controlled motors. The 
Micromegas detector has an active surface of 50x50 mm 2 and is enclosed 
in a vessel, equipped with a quartz entrance window. It is composed of 
a 1.6-mm drift gap and a 160-/Ltm amplification gap, separated by a 400- 
lpi electroformed nickel mesh. The drift electrode is composed of a 0.5-nm 
thick Ni-Cr layer deposited on a quartz lamina. When the laser is focused 
on this electrode, the high energy density produces electrons by non- linear 



photoelectric emission processes 2l|. The detector is operated with a Ne 
95 % iC4Hio 5 % mixture at atmospheric pressure with gains varying from 
10 to 10 5 . The drift field is fixed to 1 kV/cm. A 55 Fe source is placed 
inside the vessel for absolute gain calibration. Different sets of pads and 
strips of various sizes are implemented on the anode PCB. In the present 
experiment, we use a set of nine 4x4 mm 2 pads. The Micromegas mesh is 
connected to a voltage amplifier of gain 100 made with two NXP 5205A wide 
band amplifier integrated circuits [22|. The anode pads readout is achieved 



by a GASSIPLEX chip 23J], which is a 16-channel multiplexed low noise 
preamplifier shaper. It works at a gain of 2.2 mV/fC, a peaking time of 
1.2 /is and has a maximum dynamic range of 500 fC. The electronic noise 
is of 2 000 electrons RMS, which makes it possible to work at gains greater 
than 10 4 for a low number of primary electrons. Digitisation (12-bit ADC) 



zero suppression and pedestal subtraction are performed using CAEN 24] 



modules (V550A and V551). The data acquisition system is triggered either 
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by the mesh signal in runs with the Fe source, or by a Photonis XP2282B 
photomultiplier 25f] connected to the laser through an optical fibre. 



3. Gain calibration 

In principle the gain of the Micromegas detector can be deduced from 
the charge measurements in the single electron working mode. However, the 
standard calibration based on the 5.9 keV X-rays of 55 Fe has been performed 
in order to crosscheck the results in the same gain domain. Unfortunately, 
the limited dynamic range of the GASSIPLEX circuit associated with the 
relatively high number of electrons given by the source does not make it pos- 
sible to cover a wide gain range in the zone of interest. A method combining 

Micromegas 160nm -Ne 95% iC 4 H 10 5% ■ E d =1 kV/cm 
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Figure 2: Gain of the Micromegas detector in Ne 95% iC4Hio 5% as a function of the 
mesh voltage. Three sets of data are represented: 55 Fe source (cross), 10±1 triangles), 
and 20±2 (squares) electrons produced by focusing the laser on the drift electrode. 

the source and the laser has been applied and the results are presented on 
Fig. [2 The charge on the central 4x4-mm 2 pad is measured at a mesh 
voltage of 440 V with the 55 Fe source. It gives an absolute gain of 1.4 10 
with a number of primary electrons of 167. This number is the one obtained 
for a 5.9 keV X ray and for an ion-electron pair production mean energy of 
35.3 eV in the Ne 95% iC 4 Hi 5% mixture (W Ne = 36 eV and W iC4 u 10 = 
23 eV being the mean energy for Neon and Isobutane (26|). At the same 
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voltage, the laser is focused on the drift electrode in front of the central pad. 
From the gain value obtained with the source, the number of primary elec- 
trons can be deduced. Then it can be adjusted at a much lower level than 
the one provided by the source. Measurements were performed with two 
values of primary electrons (10±1 and 20±2) and gave compatible results. 
The mesh voltage has been varied from 430 V to 490 V, i.e gains from 10 
to 7.10 4 . 



4. Single-electron mode 

The method applied to reach the single electron mode is similar to the 
one used for photomultipliers 27]]. The pulsed laser light intensity is at- 



tenuated to a regime in which the probability of producing one electron is 
much greater than the probability of giving more than one. In the present 
setup, the laser light intensity is attenuated by a factor of 2 000. In those 
conditions, only 7 % of the events are considered as non-zero events (charge 
greater than 10 4 electrons, which is 5 times the RMS noise) and the proba- 
bility of obtaining more than one electron can be considered as low (0.5 %). 
The charge spectra measured on the Micromegas central pad at five differ- 
ent biases (470 to 510 V) are shown in Fig. 3. Experimental and theoretical 
studies in proportional gas counters [3] showed that the gain fluctuations 
can be described by two models depending on the electric field to pressure 
ratio R, the gas characteristics and the detector geometries. For low R, the 
probability of ionisation by an electron collision within an avalanche can 
be assumed to be independent of its past history: the SER is exponential 
according to the Furry law. For larger R values, measured SER showed pro- 
nounced maxima. This departure from the Furry distribution means that 
the ionising collision probability cannot be supposed independent of the path 
travelled from the previous ionisation: the SER follows a Polya distribution. 
Concerning MPGDs, the SER has been measured and fitted by Furry laws 
at moderate 28] and high [29! 30] gains with stacks of GEMs. It has been 



fitted by a Polya for a triple GEM working at a gain of 4.1 10 5 [29J]. The 
SER has also been measured with a Micromegas at gains greater than 10 6 
and follows a Polya distribution 3l|. Each measured SER clearly shows a 
maximum and can be fitted by a Polya distribution: 
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(a) SER Ne 95% iC H 10 5% - V Mesh =470V 



(b) SER Ne 95% iC H 10 5% - V Mesh =480V 




Polya Fit results 
X 2 /ndf: 186.7/89 
C = 144.6+/- 3.2 
Q = 29.4 +/- 0.3 
8= 2.3 +/- 0.1 
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(c) SER Ne 95% iCH 10 5% - V M , !h =490V 



Polya Fit results 
X 2 /ndf: 304.3 /156 
C = 68.1 +/- 1.3 
Q = 52.7 +/- 0.5 
b 2.3 +/- 0.1 




50 100 150 200 250 300 350 400 

Q(ch) 



(e) SER Ne 95% iC H 10 5% - V Me!h =510V 



Polya Fit results 
X 2 /ndf: 255.9 /230 
C n = 11.8+/- 0.4 
Q = 99.3 +/- 2.4 
1.7+/- 0.1 





Polya Fit results 
X 2 /ndf: 244.3 /119 
C = 104.4+/- 2.0 
Q = 39.8 +/- 0.4 
e- 2.2 +/- 0.1 



50 100 150 200 250 300 350 400 

Q(ch) 



(d) SER Ne 95% iC H 10 5% - V Mesh =500V 



Polya Fit results 
X 2 /ndf: 339.2/215 
C„= 52.3 +/-1.0 
Q„ = 71.8+/-0.7 
1.9+/- 0.1 
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Figure 3: Single-electron response (SER) of the Micromegas detector measured on one 
pad in Ne 95% iCziHio 5% . Each spectrum is fitted with a Polya function (plain line). 
The dotted line is the extrapolation to zero nf the fitted function. 



where Cq is a constant, Qq the mean charge and 9 the parameter of the 
Polya. From the mean charge Qq , a gain value is deduced for each mesh 



Mesh voltage (V) 


Gain (calibration) 


Gain (single electron) 


470 


3.7 10 4 ±0.3 


3.6 10 4 


480 


5.0 10 4 ±0.5 


4.6 10 4 


490 


7.0 10 4 ±0.7 


6.0 10 4 



Table 1: Comparisons of the gains measured by the calibration (see section 3) and in 
single-electron mode. 



voltage and for the three lower ones (470, 480 and 490 V) they are compared 
to the calibration values obtained in section 3 (see table CQ . There is a good 
agreement between the two calibrations based on the 167 electrons of the 
Fe source or deduced from the SER. Moreover, at those gains (3.7 10 4 , 5.0 
10 4 and 7.0 10 4 ), the 9 parameter does not depend on the gain: 9 = 2.2±0.2. 
Even if no direct comparison can be made (different detector and different 
gas), it is encouraging to have a result similar (9 = 2) to the one obtained 
with a triple-GEM operated in Ar±CH 4 1.3 % at a gain of 4.1 10 5 [29]. The 
measurements made at the two larger bias voltages (500 and 510 V) give 
gains which are respectively less by 10 % and 15 % outside the error bars than 
the ones expected from the fit and extrapolation of the calibration results of 
section [3) For those two gains, the parameter 9 has smaller values: 9 = 1.9 
and 1.7 for 500 V and 510 V mesh bias respectively. The discrepancies could 
be explained by non linearity in the gas multiplication process particularly 
when a low proportion of quencher is used [32]. Anyway, it has been decided 
to only retain in further studies only the three SER results for which the 
agreement with the gain calibrations is satisfactory. 

In detectors, an important feature is the gain fluctuations which can be 
expressed by the relative gain variance / of the SER. If the SER is fitted by 
a Polya distribution, the relative gain variance is simply deduced from the 



parameter 9 by 33]: 



f = TTe (2) 

From our measurements (biases at 470, 480 and 490 V) and the associated 
Polya fit, we have the same value for 9 = 2.2 ± 0.02 and thus the same for 
/ = 0.31 ±0.02. The gain variance measured in our operative conditions is 
two times lower than the one usually obtained in MWPC (/=0.7) which is 
sometimes taken as the parameter for the gain fluctuations in some MPGD's 



models [34j |. In 3JJ], the SER has been measured in a Micromegas of 100 /*m 
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amplification gap with He-Isobutane 4% and 10% gas mixtures for gains 
above 10 6 . The spectra fit well with a Polya distribution, but the gain 
variance is double of that of our present result. However our value of / is in 
good agreement with the one (/=0.33) obtained with a triple GEM operated 
in Ar+CH4 1.3% at a gain of 4.1 10 5 [29]. Let us consider our two results 
rejected for calibration considerations. The detector was working at higher 
gains (biased at 500 and 510 V) and the Polya fit gave lower values for 9 (1.9 
and 1.7) and thus larger values for / (0.34 and 0.37 respectively). It clearly 
indicates that avalanche fluctuations depend on the gain and increase at 
high gain values. That was already observed in parallel-plate proportional 
counters 35]- If the SER has to be used for energy resolution purposes, it 
must be measured at gains close to the one at which the detector is operated 
in the experiment. 

5. Determination of the energy resolution 

The energy resolution is given by the charge resolution. The relative 
charge resolution y/vai (Q)/Q is related to the number of primary electrons 
n e by: 

var(Q) _ f var(n e ) (^i\ 2 1 ,„\ 
(E(Q))2 E(n e ) + E2(n e ) \G ) E^n e ) [ > 

where / is the relative gain variance, G the detector gain and 07 the stan- 
dard deviation of the noise at the electronics input. This formula and its 
conditions of application are discussed in Appendix [Al In the current con- 
figuration, the electronic noise is 2 000 electrons RMS whereas the gain 
is always greater than 10 4 . Since / = 0.31 and (<7i/G) 2 < 0.04, the term 
((7f/(jrE(n e )) 2 is negligible compared to //E(n e ), even for low number of 
electrons. 

When the electrons are generated only by ionisation in the gas, as is the 
case with the X-rays produced by a 55 Fe source, the variance of n e can be 
expressed by the following formula: 

var(n e ) F 



E2(n e ) E(ne) (4) 

where F is a constant called the Fano factor [g, 0]. It depends on the 
interactions of the particle with the gas under consideration. 

When the primary electrons are generated by laser photons interacting 
with a metal, the variations of n e depend on the laser energy fluctuations 
and on the statistics of photon absorption: 
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var(n e ) _ ^ 2 F /as 

E2(n e ) " + E(n e ) W 

where a is a parameter related to the energy resolution of the laser, and. Fias 
is equivalent to a Fano factor associated with the electron emission from the 
photons. The details of the calculation and the conditions of application are 
also given in Appendix [XI 

These expressions of the variance of n e have been used to characterise 
our Micromegas detector operating with Ne 95 % iC4H 10 5 % gas mixture. 

5.1. Measurement with the laser 

The expression of the charge resolution measured with the laser is ob- 
tained by applying ([5j) to (J3J) : 



Q J E(Q) V E ( n 
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« 80 
c 

3 

o 

O 70 



10 



Pad 

Entries 2000 
Mean 0.1985 
RMS 0.03525 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

Q(pC) 



Figure 4: Charge induced on one pad of the Micromegas detector by 33±3 electrons 
produced by focusing the laser on the drift electrode. The gas mixture is Ne 95% iC4Hio 
5%, and the mesh voltage 470 V. 
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Micromegas 160n -Ne 95% iC^H,,, 5% - Energy resolution 




t 0.1 0.2 0.3 0.4 0.5 0.6 

Figure 5: Energy resolution as a function of the number of primary electrons produced by 
focusing the laser on the drift electrode and varying the light intensity. The data points 
are fitted following a quadratic function (plain line). The dotted line is an extrapolation 
of the fit. The point corresponding to the Fe source is added on the graph. 

An example of charge spectrum obtained with 33±3 primary electrons 
is shown in Fig. [H The formula ([6]) was used to fit the charge resolution 
versus the number of injected electrons, as shown in figure [5j The actual 
fit function, which is shown in the figure inset, uses the parameters a and 
b, which stands for \/f + Fi as . A very good agreement is observed. A 
value of 0.71 ± 0.10 for Fi as is deduced from the fit and the value of / 
determined from the single electron response measured previously. The a 
value (3 ± 1 %) is consistent with the laser fluctuations measured with 
the control PMT (2 %), as well as the value given by the manufacturer (< 
4 %). Its contribution is negligible on the main part of the studied range. 
Once the laser Fano factor has been determined, the dependence of the 
energy resolution from the mesh voltage was investigated. Table [2] sums up 
the results. It shows that the energy resolution is independent of the gain 
in the voltage range of operation. This result is consistent with the SER 
measurements presented in section [H 
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Number of primary electrons n e 


Mesh voltage (V) 


Gain 


Charge resolution 






^ 7 1 D 4 






4oU 


c n iM 
O.U 1U 


n ap. 

U.4D 




490 


7.0 10 4 


0.44 


lOrbl 


440 


1.4 10 4 


0.32 




470 


3.7 10 4 


0.35 




480 


5.0 10 4 


0.34 




490 


7.0 10 4 


0.32 


20±2 


440 


1.4 10 4 


0.23 




470 


3.7 10 4 


0.23 




480 


5.0 10 4 


0.22 


33±3 


440 


1.4 10 4 


0.17 




470 


3.7 10 4 


0.19 



Table 2: Charge resolution for several number of primary electrons at different gains 



5.2. Measurement with a 55 Fe source 

In the case of a particle impinging the gas, the average number of primary 
electrons is defined by: 

EK) = L 

where W is the mean ionisation energy and £ the energy of the incident 
particle (see section [3] for the determination of their values). Then, the 
expression of the energy resolution measured with a 55 Fe source is obtained 
by applying g]) to ©: 



Q J E(Q) V w ' £ 

The charge spectrum measured on the central 4x4- mm 2 pad at a mesh volt- 
age of 440 V is shown in Fig. [6j The events were selected by anticoincidence 
with the neighbouring pads. The value of the relative gain variance / for 
a mesh bias of 440 V is the same than at 470, 480 and 490 V, as shown 
in section 15.11 From the measured charge resolution, the value of the Fano 
factor deduced is 0.37. Due to the low counting rate on the central pad, 
the data had to be acquired over almost 10 hours, and the gas gain could 
change during the experiment. If such a variation occurs, the actual value of 
/ is increased and the measured spectrum is wider than expected. Thus the 
value given for the Fano factor can only be interpreted as a maximum value: 
F < 0.37. This result is a first experimental estimation of this parameter 
for a binary mixture consisting of Neon and a molecular gas. 
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Source 55 Fe- E d =1kV/cm - V Mesh =440V 
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Fit results 
a„ + a,Q + A exp(-[(Q-Q o )/af 

X 2 /ndf : 78.2 /56 
a = 57.33 +/- 3.35 pC 
a, = -128.17+/- 8.79 
Q = 0.374 +/- 0.001 pC 
a = 0.024 +/- 0.001 pC 
a Q IQ (%): 6.4 ± 0.4 
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Figure 6: Charge induced on one pad by 5.9 keV X-ray of a Fe source. The gas mixture 
is Ne 95% iCiHio 5% and the mesh voltage of 440 V. The spectrum is fitted following 
a function (plain line), which superposes a linear function (background) and a Gaussian 
(source peak). 
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6. Conclusion 

We succeeded in separately measuring the contribution to the energy- 
resolution of the gain and the ionisation fluctuations of a Micromegas de- 
tector. The method relies on a test-bench made of a 337-nm wavelength 
laser and an optical focusing system. The detector was operated in Ne 95 % 
iC4Hio 5 %. The single-electron response (SER) was measured for gains 3.7 
10 4 , 5.0 10 4 and 7.0 10 4 . The experimental data were fitted with a Polya 
distribution. The obtained Polya parameter ( = 2.2 ± 0.2) allowed us to 
deduce the relative gain variance (/ = 0.31± 0.02) which did not depend 
on the gain in the range of operation. The laser system enabled us to study 
the energy resolution as a function of the number of primary electrons n e . 
A Fano factor could be deduced in the case of the laser-metal interaction: 
Fias = 0.71 ± 0.10, taking the value of / obtained from the SER. A first ex- 
perimental determination of the maximum value of the Fano factor for the 
Ne 95% iC4Hio 5% mixture at atmospheric pressure with 5.9 keV X-rays 
was obtained: F < 0.37 . 

The experimental method presented here could be widely used for vari- 
ous MPGD's geometries and gas mixtures. It could give experimental mea- 
surements of the gain variance and the Fano factor, which are fundamental 
parameters of these detectors. The detector used in the present work showed 
good performances in terms of energy resolution at atmospheric pressure, 
thanks to its low gain fluctuations. Further studies of its properties at dif- 
ferent pressures and gas mixtures are planned in the future. 
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A. Charge resolution 

The aim of this appendix is to demonstrate the formula ([3]) used in 
section [5j 

It is assumed that there is no attachment of primary electrons during the 
drift, that the multiplication in the gas has linear behavior, that all multi- 
plications are independent from one another and have the same probability 
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law. In this case, the response to n e electrons is the sum of n e single electron 
responses. Then, the number Q of electrons detected at the anode is: 



Q = N noise + Y, N SER,k (A.l) 

l<k<n e 

where N$ er is the random variable of electrons collected at the anode for a 
single photoelectron (single electron response or SER) and N no i se the ran- 
dom variable corresponding to the noise. The sum is calculated over a 
random number of terms, n e . When n e = , Q is evaluated only from 
the noise. The expectation and variance of Q are calculated from equation 
(|A.1|) assuming that the noise is totally independent from the multiplication 
process: 

E(Q) = E(N noise ) + E(N SER )E(n e ) 

var(Q) = vav(N noise ) + var(NsER)E(n e ) + E 2 (iVs£; fl )var(n e ) 

The demonstration can be achieved with conditional probabilities, which is 
less demanding on the properties of the distribution than the method of 
moment generating functions. The average of the noise is usually 0. Let / 
designate the relative variance of the single electron response: 

_ var(NsER) 
1 ~ E 2 (N SER ) 
and define the following variables: 

a i = \f vai(N noise ) 
G = E(Nser) 

corresponding to the standard deviation of the noise an to the detector gain 
respectively. Then: 

var(Q) _ f var(n e ) (Qi\ 2 1 , . ^ 

(E(Q))2 E(n e ) + E2( ne ) + lGi E^) { ' } 

In the case of a radioactive source, var(n e )/E 2 (n e ) can be expressed with 
the introduction of the Fano factor 0, 0] : 

var(rip) F 



E 2 (n e ) E(n e ) 

In the case of the laser, var(n e )/E 2 (n e ) depends on the laser light fluctua- 
tions and on the electron emission statistics, which can be expressed by: 
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var(n e ) 



Flas 

E(n e ) 



E 2 (n e ) 



The second term i^ as /E(n e ) expresses the contribution of the fluctuations 
in the electron emission process and has been modelled as a deviation from 
a standard Poisson law by the introduction of the parameter F[ as , which is 
analogous to a Fano factor. The a parameter is related to the laser light 
fluctuations. As discussed in Ref. [211 ] . the electron emission follows a power 
law with an exponent p of around 2 . For the small variations of the laser 
light q (around 2 % in our setup), a linear expansion of the power law leads 
to multiply a by p and thus a = pa ~ 4 %. Equation (|A.2p then becomes: 
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